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APPLICATION  OF  2 FRACTIONAL  FACTORIALS  IN  SCREENING 
OF  VARIABLES  AFFECTING  THE  PERFORMANCE  OF  DRY  PROCESS 

ZINC  BATTERY  ELECTRODES 

Nicholas  T.  Wilburn 

U.  S.  Army  Electronics  Research  and  Development  Laboratory 
Fort  Monmouth,  New  Jersey 

Among  its  research  and  development  activities  on  zinc-silver  oxide 
batteries  for  special  applications,  USAELRDL  is  investigating  the  prep- 
aration of  zinc  electrodes  by  dry  processes.  These  involve  the  appli- 
cation of  dendritic  zinc  powders  under  pressure  to  the  grids.  The  inter- 
locking properties  of  the  dendritic  zinc  particles  make  it  possible  to  form 
the  electrodes  with  moderate  pressures  such  that  the  porosity  and  re- 
lated high  surface  area  of  the  electrode  is  not  destroyed.  It  is  expected 
that  dry  process  zinc  electrodes  will  have  many  advantages  over  conven- 
tional electrodepo sited  sponge  zinc  electrodes,  including  higher  discharge 
efficiency,  greater  uniformity  of  performance  and  better  adaptability  to 
mechanized  production  with  resultant  economics. 

Due  to  the  large  number  of  variables  affecting  the  discharge  perform- 
ance of  the  electrodes,  it  was  decided  to  design  and  conduct  a fractional 
factorial  experiment  to  isolate  the  significant  variables.  These  variables, 
and  any  controlling  interactions  between  them,  could  then  be  studied  fur- 
ther to  arrive  at  the  optimum  conditions  for  the  production  of  electrodes 
of  maximum  discharge  efficiency.  The  fractional  factorial  experiment 
was  thus  intended  for  the  preliminary  screening  of  all  major  variables 
acting  simultaneously.  As  such,  it  was  recognized  that  it  is  the  most 
efficient  and  economical  process  known  for  accomplishing  this,  in  addit- 
ion to  providing  valuable  data  on  interactions  between  the  variables 
which  cannot  be  obtained  from  the  more  widely  used  one  or  two  at  a time 
variable  investigations. 

There  were  two  categories  of  variables  in  the  electrode  investigation, 
those  related  to  the  electrolytic  formation  of  the  dendritic  zinc  powders 
and  those  related  to  the  electrode  preparation  itself.  Although  several 
more  variables  were  considered,  it  was  decided  to  limit  the  number  of 
variables  to  eight,  keeping  all  other  factors  constant.  The  eight  selected 
variables  are  shown  in  Figure  1.  High  and  low  levels  as  shown  were 
assigned  to  each  variable.  Although  considerable  thought  was  given  to 
determination  of  the  levels,  .it  is  seen  in  retrospect  that  wider  ranges 
might  have  been  assigned  in  some  instances.  Based  on  available  literature, 
these  were,  however,  considered  sufficiently  wide  ranges. 

Figures  can  be  found  at  the  end  of  this  article. 
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Variable  G,  the  electrolyte  temperature  during  the  plating  operation, 
turned  out  to  be  impossible  to  control  with  the  plating  equipment  which 
was  prepared  for  the  experiment  and  with  the  selected  plating  current 
densities.  The  experiment  thus  became  one  involving  seven  variables, 
each  at  two  levels.  Rather  than  to  eliminate  G,  however,  it  was  felt 
that  this  paper  would  serve  a broader  purpose  if  it  gave  the  fractional 
factorial  procedures  for  the  study  of  eight  as  well  as  seven  variables. 
Variable  G should,  therefore,  be  considered  both  in  and  out  of  the 
figures  and  analysis. 

Having  established  the  variables  and  the  high  and  low  levels,  the 
fractional  replicate  design  was  established,  as  shown  in  Figure  2.  This 
is  the  first  design,  a second  having  been  run  later  for  reasons  to  be  dis- 
cussed. The  design  involves  eight  variables  (or  seven)  each  at  two  levels 

g 

to  be  studied  with  a total  of  sixteen  different  electrodes.  A full  2 fac- 
torial experiment,  eight  variables  each  at  two  levels,  would  involve  2® 
or  256  trials  (128  trials  for  a 2 factorial).  Therefore,  this  design 

Q 

represents  a one  sixteenth  fraction  of  the  2 factorial  (one  eighth  of  the 
2*^),  expressed  as  a 2®"^  design  (or  27-3  for  seven  factors).  The  de- 
sign is  based  on  extension  of  a basic  2^  or  16  trials.  Thus  the  first 
four  variables  are  arranged  in  standard  order  for  the  2^  factorial.  The 
other  variables  are  then  introduced  by  making  a basic  assumption  that 
three  factor  interactions  between  the  first  four  variables  are  negligible. 
Therefore  E is  introduced  by  equating  it  to  the  interaction  between 
variables  A,  B,  and  C.  Regarding  the  high  level  as  plus  and  the  low 
level  as  minus,  the  level  of  E for  the  first  box  is  - x-  x-  = -, 

+ X - X - = + for  the  second  box,  etc.  Similarly,  variable  F is  intro- 
duced to  equating  it  to  the  BCD  interaction,  variable  G to  the  ABD 
interaction  and  variable  H to  the  A CD  interaction.  Normally  in  a seven 
variable  design,  the  seventh  variable  would  be  equated  to  ABD  instead  of 
A CD  as  was  done  here  since  variable  G was  dropped  out.  However, 
this  does  not  affect  the  experiment  in  any  way. 

Having  thus  established  the  fractional  design,  the  sixteen  electrodes 
were  prepared  in  accordance  with  the  high  and  low  level  criteria  for  each 
variable.  The  electrodes  were  then  tested  one  at  a time  under  standard 
and  carefully  controlled  conditions  to  give  the  sixteen  yields  or  responses 
as  shown.  The  response  is  the  discharge  efficiency  of  the  electrode  ex- 
pressed in  percent  as  the  ratio  of  the  output  capacity  (at  a discharge 
current  density  of  1.  1 amp/ square  inch  to  a 0.  3 volt  change  for  the  elec- 
trode) to  the  theoretical  capacity  of  the  zinc  active  material.  For  com- 
parison the  conventional  sponge  zinc  electrodes  give  an  average  efficiency 
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of  about  20%  under  comparable  discharge  conditions.  The  average  of  the 
sixteen  responses  is  31.  2%.  Before  proceeding  with  the  analysis  of  the 
responses,  it  should  be  noted  how  the  design  can  be  used  to  relate  the 
sixteen  mean  effects  obtained  from  the  response  analysis  to  the  variables 
and  their  two  factor  interactions.  Consider  the  high  level  boxes  in  the 
first  four  columns.  The  first  mean  effect,  no  high  levels,  will  be  for 
twice  the  average.  The  second  will  be  for  variable  A,  the  third  for  B , 
the  fourth  for  the  AB  interaction,  the  fifth  for  C and  so  on.  The  eighth 
(A,  B and  C at  the  high  level)  will  be  for  variable  E which  was  origin- 
ally equated  to  ABC.  Similarly  the  twelfth  for  G,  the  fourteenth  for  H 
and  the  fifteenth  for  F.  Since  this  is  a one  sixteenth  fraction,  each  of 
these  principal  effects  will  be  confused  with  fifteen  other  effects.  How- 
ever, these  for  each  of  the  eight  variables  will  be  three  factor  or  higher 
order  interactions  which  are  considered  negligible,  the  basis  on  which 
the  design  was  established.  Each  two  fe.ctor  interaction,  AB  for  example, 
will  be  confused  with  fifteen  other  effects  of  which  three  are  other  two 
factor  interactions.  The  sixteenth  row  of  the  First  Design  will  represent 
such  a combination  of  four  2 factor  interactions. 

The  analysis  of  the  sixteen  responses  is  shown  in  Figure  3.  The 
technique  used  here  is  the  Yates'  Algorithm  which  is  a rapid  method  for 
obtaining  the  same  mean  effects  that  would  be  obtained  from  a formal 
and  lengthy  analysis  of  variance.  The  Yates'  Algorithm  is  applicable 
to  any  factorial  experiment.  Its  advantages  become  more  apparent  the 
larger  the  experiment. 

The  mechanics  of  the  Yates'  calculations  are  very  simple.  The  first 
figure  in  Column  (1)  is  the  sum  of  responses  1 and  2;  the  second,  the 
sum  of  responses  3 and  4,  etc.  The  ninth  figure  is  the  sum  of  responses 
1 and  2 with  the  sign  of  response  1 reversed.  Column  (2)  is  derived  from 
Column  (1)  in  the  same  manner.  Additional  columns  are  introduced  until 
a column  is  completed  the  first  figure  of  which  is  equal  to  the  sum  of  the 
responses.  The  arithmetic  in  each  column  is  checked  before  proceeding 
to  the  next  column.  The  sum  of  Column  (1)  is  equal  to  twice  the  sum  of 
the  even  numbered  responses;  the  sum  of  Column  (2)  is  equal  to  four 
times  the  sum  of  every  second  even  numbered  response;  the  sum  of 
Colunon  (3)  is  equal  to  eight  times  the  sum  of  every  fourth  even  numbered 
response;  and  the  sum  of  Column  (4)  is  equal  to  sixteen  times  the  sixteenth 
response.  The  mean  effects  are  obtained  by  dividing  each  figure  of  Column 
(4)  by  eight.  The  sum  of  the  mean  effects  is  also  checked.  The  62.  363 
figure,  twice  the  average,  is  not  used  in  the  subsequent  analysis.  The 
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effects  measured  by  the  mean  effects  are  given  in  the  last  column.  The 
effects  A»  B,  AB,  C,  AC,  etc.  are  those  as  previously  read  off  the  first 
design  chart.  The  three  factor  and  higher  order  interactions  which  they  are 
confused  with  are  of  no  importance  since  the  design  is  based  on  their 
being  assumed  negligible.  The  other  two  factor  interactions  are  impor- 
tant and  must  be  known.  They  are  obtained  from  the  effects  A,  B,  AB, 
etc.  and  the  defining  contrasts  of  the  design.  These  defining  contrasts 
are  obtained  from  the  equating  that  had  been  done  in  establishing  the 
fractional  design;  E = ABC,  F=  BCD,  G = ABD  and  H = ACD.  Sixteen 
defining  contrasts  are  required  for  the  eight  variable  design.  The  first 
defining  contrast  is  always  I,  the  next  four  are  ABCE,  BCDF,  ABDG 
and  ACDH.  The  remaining  eleven  are  found  by  exhaustively  multiplying 
these  contrasts  using  the  rule  that  like  factors  are  cancelled  out. 


(1)  I 

(2)  ABC  X E = ABCE 

(3)  BCD  X F = BCDF 

(4)  ABD  X G = ABDG 

(5)  ACD  X H = ACDH 

(6)  A^CE  X ^<XDF  = ADEF 

(7)  A^CE  xAiiDG  = CDEG 

(8)  AB^E  xA<?DH  = BDEH 

(9)  SCUF  X A1^0G  = ACFG 

(10)  BC0F  X AC0H  = ABFH 

(11)  AB0G  X ACJ&H  = BCGH 

(12)  A^<?E  X )^C0F  X AB0G  = BEFG 

(13)  AWE  X IJC0F  X AO^H  = CEFH 

(14)  A1S(7E  X A00G  X ACi^H  = AEGH 

(15)  1^(^0F  X A^l^G  xACDH  = DFGH 

(16)  A0CE  X ]^C0F  X ABlAG  X ACDH  = ABCDEFGH 


(These  sixteen  defining  contrasts  may  also  be  obtained  by  simply 
reading  off  the  low  level  boxes  for  each  of  the  sixteen  trials  in  the  First 
Design  chart.  However  this  procedure  will  not  apply  in  all  cases,  e.  g. 
in  the  2^“^  it  would  give  sixteen  defining  contrasts  when  only  eight  are 
required.  ) 

The  principal  effects  are  found  by  multiplying  the  first  effect,  A,  B, 
AB,  etc.  , by  the  sixteen  defining  contrasts.  The  two  factor  interactions 
confused  with  AB  are  found  for  example  to  be  CE,  DG  and  FH. 
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AB  X I = AB 
Ati  xjSC^CE  =CE 
At  xiJf^DG  = DG 
At  xA^FH  = FH 


The  other  twelve  effects  confused  with  AB,  all  three  factor  or 
higher  order  interactions,  could  be  found  with  the  other  twelve  defining 
contrasts  if  desired.  The  total  256  effects  found  by  multiplying  each 
starting  effect  by  the  sixteen  defining  contrasts  would  give  the  full  256 

effects,  from  I to  the  interaction  ABCDEFGH,  which  would  be  obtained 

8 

in  conducting  a full  2 ^ctorial  experiment. 

7-3 

The  defining  contrasts  for  the  2 experiment  are  I,  ABCE, 

BCDF,  ACDH,  AEDF,  BDEH,  ABFH,  and  CEFH.  The  effects  measured 
are  found  in  the  same  way  as  for  the  eight  variable  design.  The  only 
difference  is  that,  in  this  case,  all  effects  involving  G drop  out.  The 
twelfth  set  of  effects  measured,  identified  by  the  asterisk,  then  becomes 
eight  interactions,  all  three  factor  or  higher  order. 

The  fifteen  mean  effects  with  their  identifying  effects  are  then  ordered 
by  arranging  them  in  order  of  magnitude  without  regard  to  sign.  The 
thus  ordered  set  of  mean  effects  is  then  arranged  in  a half -normal  plot 
as  shown  in  Figure  4 to  interpret  the  relative  significance  of  the  effects. 
The  ordinate  is  the  order  number  of  the  fifteen  effects  from  smallest  to 
largest.  The  abscissa  gives  the  mean  effect  magnitudes.  The  fifteen 
points  which  are  plotted  are  identified  by  the  proper  major  effects  and 
two  factor  interactions.  The  plot  is  given  for  both  the  full  eight  variables, 
and  the  actual  seven  variable  experiment.  In  the  latter  case,  the  G 
factor  and  the  G interactions  drop  out.  The  asterisk  again  denotes  high 
order  interactions.  In  the  plot  an  error  best  straight  line  has  been  drawn 
through  the  lowest  seven  points.  High  magnitude  effects  falling  signifi- 
cantly off  the  line  are  judged  to  be  distinct  from  error  and  therefore  con- 
trolling factors  in  the  process  being  investigated.  In  a plot  o.f  this  type, 
it  is  considered  unusual,  however,  to  have  a well  defined  errbr  line 
with  as  many  as  eight  points  falling  clearly  off  it.  To  gain  insight  into 
this  unusual  behavior  as  well  as  to  gain  more  precision  in  the  estimation 
of  all  of  the  effects,  it  was  decided  to  conduct  a second  phase  of  the  experi- 
ment, another  group  of  sixteen  electrodes  differing  from  the  first.  Since 
an  interaction  between  A (the  highest  magnitude  effect)  and  B appeared 
reasonable,  and  since  F and  H were  both  high  magnitude  effects  which 
might  interact  with  each  other,  it  was  decided  to  establish  the  second 


22 


Design  of  Experiments 


design  in  such  a way  as  to  separate  the  AB  and  FH  interactions. 

Although  this  can  be  done  in  several  equally  effective  ways,  it  was  de- 
cided to  do  it  by  reversing  the  levels  for  factors  A and  D,  giving  the 
design  as  shown  in  Figure  5. 

This  design  is  identical  to  the  first  except  for  the  reversal  of  levels 
of  variables  A and  D.  Sixteen  electrodes  were  prepared  in  accordance 
with  the  indicated  variable  levels.  The  electrodes  were  then  tested  to 
give  the  listed  responses.  The  next  step  in  the  experiment  was  to  deter- 
mine the  mean  effects  generated  by  these  responses  and  to  combine  them 
with  the  mean  effects  resulting  from  the  first  design.  This  may  be  done 
in  two  ways.  The  first  is  to  combine  the  two  sets  of  sixteen  responses 
into  an  overall  group  of  thirty-two  and  then  to  conduct  a Yates'  computat- 
ion to  arrive  at  the  thirty-two  mean  effects.  The  second  way  is  to  conduct 
a Yates'  computation  on  the  second  design  responses  and  then  combine 
the  mean  effects  with  those  of  the  first  design.  The  first  method  is  less 
time  consuming  and  therefore  preferable.  The  second  method  gives  a 
clearer  picture  of  the  separation  of  effects  and  is  therefore  now  given. 

(The  first  method  is  given  in  Appendix  A-1) 

Figure  6 shows  the  Yates'  computations  on  the  responses  from  the 
second  design.  The  operations  are  identical  to  those  as  described  for 
the  first  design.  It  is  noted  that  the  signs  of  all  elements  involving  A 
and  D in  the  Effects  Measured  column  are  now  minus  since  the  levels 
of  A and  D had  been  reversed.  This  reversal  of  signs  permits  the 
separation  of  effects  as  shown  in  Figure  7. 

The  mean  effects  derived  from  the  first  design  are  given  in  the  column 
marked  X.  Those  derived  from  the  second  design  are  given  in  the  column 
marked  Y.  An  example  of  the  computation  to  separate  the  effects  is  as 
follows: 

The  second  mean  effect  in  the  X colvunn  is  for  variable  A.  The 
similar  mean  effect  in  the  Y column  is  for  minus  variable  A with  the 
difference  in  the  absolute  magnitudes  of  the  mean  effects  being  due  to 
experiihental  error.  Reversing  the  sign  of  the  column  Y mean  effect 
and  averaging  it  with  that  of  the  column  X mean  effect  will  give  the 
-4.  81  mean  effect  for  variable  A and  certain  high  order  interactions 
confused  with  it.  1/2(X+Y),  -1.  36,  then  gives  the  mean  effect  for  the 
remaining  high  order  interactions  originally  confused  with  A.  Similar 
calculations  are  performed  to  separate  all  of  the  other  effects,  those 
containing  A or  D from  each  of  the  others.  The  thirty-one  statistics 
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thus  obtained,  not  including  the  61.22  (twice  the  average)  figure,  are 
then  arranged  in  order  of  magnitude  without  respect  to  sign  and  plotted 
in  a thirty-one  factor  half -normal  plot  as  shown  in  Figure  8. 

It  is  seen  that  the  error  best  straight  line  is  now  established  by 
twenty-three  of  the  thirty-one  points  reflecting  the  greater  precision 
achieved  by  doubling  the  experiment  and  thereby  reducing  the  variance 
of  the  estimated  effects  by  one -half.  Of  the  eight  points  which  are 
clearly  off  the  line,  two  of  them,  denoted  by  asterisks, are  combinations 
of  high  order  interactions.  Their  relative  significance  cannot  be  inter- 
preted within  the  limits  of  the  experiment.  As  will  be  seen,  equating 
them  to  zero  will  not  affect  the  results.  The  BF  interaction  is  not 
far  off  the  line  and  its  significance 'may  be  questioned.  Variables  E 
and  A are  both  clearly  controlling  factors  in  the  efficiency  of  the  zinc 
electrodes.  The  signs  of  their  mean  effects  are  both  minus,  indicating 
that  higher  efficiencies  can  be  obtained  at  the  lower  levels  of  the  ranges 
studied,  in  other  words  at  the  lower  pressing  temperature,  80°  F, 
and  with  the  smaller  weight  of  zinc  per  plate.  The  interpretation  of 
variables  H,  the  formation  current  density,  and  F,  the  presence  or 
absence  of  zinc  oxide  in  the  formation  electrolyte,  is  more  complex. 
This  results  from  the  probable  significance  of  FH,  the  interaction  be- 
tween them  (CE  is  very  unlikely  to  be  significant  due  to  the  low  magni- 
tude of  C).  In  general  when  an  interaction  is  large,  as  in  this  case,  the 
corresponding  mean  effects  cease  to  have  much  meaning.  The  effect  of 
F is  clearly  dependent  upon  the  level  of  H and  vice  versa.  The  three 
effects  F,  H and  their  interaction  FH  may  best  be  interpreted  as  a 
single  highly  significant  effect.  Further  experimental  work  at  inter- 
mediate levels. for  the  two  variables  is  definitely  indicated. 

The  final  step  in  the  analysis  was  to  determine  if  the  conclusion  was 
correct  that  only  the  effects  E,  A and  the  combined  effects  of  H,  F and 
FH  were  significant,  i.  e.  , distinct  from  experimental  error.  Part  of 
the  purpose  of  this  final  step  was  to  determine  if  the  entire  experiment 
was  valid,  in  other  words,  that  there  were  no  large  errors  made  in  the 
actual  responses  which  could  have  seriously  altered  the  mean  effects. 
The  procedure  used  was  to  determine  the  standard  error  of  the  individ- 
ual observed  responses  by  analyzing  the  thirty-one  mean  effects.  A 
second  standard  error,  for  the  differences  between  observed  and  pre- 
dicted responses,  was  then  obtained  with  the  predicted  responses  based 
on  the  assumption  that  all  mean  effects  other  than  those  for  E,A,  H, 

FH  and  F were  indeed  zero.  If  the  two  standard  errors  would  then  be 
equivalent,  both  the  total  experiment  and  the  conclusions  derived  from 
it  would  be  proved  valid. 
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Since  the  error  straight  line  on  the  final  half-normal  plot  was  es- 
tablished by  the  twenty-three  lowest  magnitude  points,  a half-normal  plot 
for  these  points  was  prepared.  The  standard  error  for  the  individual 
observed  responses  as  derived  from  this  plot  was  3.  0.  (See  A-5).  All 
mean  effects  other  than  the  twice  the  average  effect,  E,  A,  H,  FH  and  F 
were  then  equated  to  zero  and  a reverse  Yates'  computation  was  con- 
ducted to  obtain  the  thirty-two  predicated  responses  (See  A -2).  These 
responses  were  compared  with  the  observed  responses  and  a list  of  the 
thirty-two  differences  between  the  observed  and  predicted  responses  was 
prepared  (See  A -3).  The  magnitudes  of  the  differenced  were  plotted  on 
normal  probability  paper.  A standard  error  was  obtained  from  this 
plot  for  the  difference  between  individual  observed  responses  and  individ- 
ual predicted  responses.  The  value  of  this  standard  error  was  3.  2 (See 
A -4).  This  was  in  excellent  agreement  with  the  standard  error  of  the 
individual  observed  responses,  thus  proving  the  validity  of  the  experiment 
and  the  conclusions  derived  from  it. 

In  conclusion,  a fractional  factorial  designed  experiment,  involving 
two  2^"^  fractional  designs,  has  been  conducted  to  determine  the  signi- 
ficance of  seven  major  variables,  and  two  factor  interactions  between  them, 
on  the  discharge  efficiency  of  the  dry  process  zinc  battery  electrodes. 

A total  of  thirty-two  electrodes  was  prepared  and  tested.  Analysis  of 
their  responses  has  indicated  the  controlling  influence  of  two  of  the  vari- 
ables, pressing  temperature  and  the  amount  of  zinc  per  plate,  and  of  the 
interaction  between  two  other  variables  relating  to  the  plating  conditions 
under  which  the  zinc  material  was  prepared.  The  other  three  variables 
have  been  shown  to  be  unimportant  in  comparison,  within  the  range  of 
levels  selected.  The  experiment  has  hilfilled  its  basic  purpose,  narrowing 
down  the  range  of  variables  to  permit  extensive  investigation  of  the  truly 
important  variables  in  order  to  arrive  at  the  optimum  electrode  prepara- 
tion procedures  in  the  most  expeditious  manner. 
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Figure  5 
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SECOND  DESIGN 
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SECOND  DESIGN  ANALYSIS,  lAOES'  ALGORITHM 
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Half -Normal  Plot , Combined  Designs 
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APPENDIX 

This  appendix  contains  supplemental  data  as  follows: 

A -1.  Yates'  Algorithm  computation  for  thirty-two  responses,  combin- 
ation of  Designs  1 and  2. 

A -2.  Reverse  Yates'  Algorithm  computation  for  thirty-two  effects,  assum- 
ing  all  effects  are  zero  except  those  for  average,  E,  A,  H,  FH  and  F. 

A - 3.  Comparison  of  observed  and  predicted  responses. 

A -4.  Probability  plot  to  obtain  standard  error  of  individual  differences 
between  observed  and  predicted  responses. 

A-5.  Half-normal  plot  of  twenty-three  mean  effects  to  obtain  standard 
error  of  individual  observed  responses. 
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A-1.  Yates'  Algorithm  computation  for  thirty-two  responses,  combination  of 
Designs  1 and  2. 
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A-3.  Ccmparison  of  observed  and  predicted  responses. 
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